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Random band matrices relevant for open chaotic systems 
are introduced and studied. The scattering model based on 
such matrices may serve for the description of preequilibrium 
chaotic scattering. In the limit of a large number of open 
channels we calculate the average 5-matrix and S-matrix's 
pole distribution which are found to reduce to those of the 
full matrix (GOE) case under proper renormalization of the 
energy scale and strength of coupling to the continuum. 



1. It is now generally acknowledged that the Random 
Matrix Theory (RMT) provides an adequate tool for the 
description of statistical properties of chaotic quantum 
systems [Q- The most of works in the RMT deals 
with the Hermitian matrices which are, strictly speaking, 
appropriate for the description of closed systems. Mean- 
while, any exited state of a quantum system has a finite 
lifetime, decaying eventually into open channels. It is 
well known |^]- ||^ that the openness of a system can be 
incorporated into consideration by means of the effective 
nonhermitian Hamiltonian 
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Due to coupling to the continuum, the effective Hamil- 
tonian (|l|) acquires, apart from the intrinsic part H, the 
antihcrmitian part which consists of the sum of products 
of the transitions amplitudes between N internal (|n)) 
and M open channel (|c)) states. For the T-invariant the- 
ory these amplitudes as well as H can be chosen to be 
real, the matrix Ti. being symmetric. The eigenvalues 
£n = En — of TL are the complex energy levels of an 
unstable system, with and r„ being, respectively, the 
energy and width of the n-th level. 

Assuming the intrinsic dynamics to be chaotic, the her- 
mitian part of TL is usually supposed to belong to the 
Gaussian Orthogonal Ensemble (GOE) of random matri- 
ces. As to the coupling amplitudes, they are considered 
to be either fixed |^,^ or random . Resulting general- 
ization of the RMT on the unstable systems proved 
pO[ | to be useful for the description of statistical proper- 
ties of quantum chaotic scattering. The resonance part 
of the scattering matrix is represented (see, e.g. [B) as 
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in terms of the effective Hamiltonian (|i|) which describes 
the evolution of the unstable system formed at the in- 
termediate stage of the collision. Therefore, fluctuations 
in scattering reflect statistical properties of the complex 
energy levels of this intermediate system. 

Recently, the considerable progress has been achieved 
in the study of dynamical and statistical properties of S- 
matrix's poles and their connection to those of scattering 
- ||l^. The limit of very large number of equivalent 
channels has a special interest since it can be related to 
the semiclassical approach [|ll|,^ . The noteworthy prop- 
erty of the pole distribution was found |^,^ in this limit: 
there always exists a finite gap between the upper edge of 
the distribution of complex energies of resonances and the 
real energy axis. This gap turns out to be the important 
characteristic of local fluctuations in chaotic scattering 
§0. 

The Gaussian distribution of the matrix elements of 
the internal Hamiltonian H implies the invariance with 
respect to the choice of the intrinsic basis. However, the 
mean-field basis seems to play an exceptional role in 
many-body chaotic systems. The realistic Hamiltonian 
matrix has a banded structure in such a representation 
|l|,g. In fact, Wi gner was the first who considered band 
matrices in connection with properties of (stable) compli- 
cated systems |]l5| . He assumed the nonzero matrix ele- 
ments to be Gaussian random variables. Now, the theory 
of random band matrices (RBM) attracts a great inter- 
est and is claimed to be relevant for studying quantum 
chaos (for review see [|l^). The most of known analytical 
results for RBM has been recently obtained by Fyodorov 
and Mirlin |17 1^ who studied RBM in the context of 
the localization, having reduced the RBM problem to the 
Efetov's supersymmetric cr-model [T9| ] with the diffusion 
constant proportional to the square of the bandwidth. 

The purpose of this Letter is to extend the RBM ap- 
proach to the consideration of open chaotic systems de- 
scribed by the effective Hamiltonian (|l]). The scattering 
model based on such random matrices can be related to 
preequilibrium chaotic scattering in a sense intermediate 
between fully chaotic scattering described by the GOE 
model and multistep compound reactions In the 

"pure" GOE models all internal degrees of freedom be- 
ing uniformly involved, the intermediate unstable system 
has already attained its complete thermodynamical equi- 
librium before a decay takes place. The band structure of 
H leads to the localization of intrinsic wave functions, an 
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intermediate decaying state concluding only a part of the 
degrees of freedom. Therefore, in addition to the decay 
time, the new diffusion time scale appears in the problem 
which characterizes the internal relaxation time. In this 
Letter we calculate the average ^-matrix and S'-matrix's 
pole distribution in the limit of very large number of open 
equivalent channels. 

2. We suppose that the hermitian part of (|l|) belongs 
to the Gaussian RBM ensemble defined by 
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where the function Jnm = J{\n — m\) decreases suffi- 
ciently fast then |n r7i| > b, with 6^1 being the 
effective bandwidth. In the GOE case, when all matrix 
elements Jnm are equal to j^, 2A determines the radius 
of Wigner's semicircle low of the eigenvalue distribution 
[Q . The transition amplitudes are considered as fixed 
quantities subject to the condition 



(4) 



which turns out to be enough for excluding direct re- 
actions. The dimensionless parameter 7 characterizes 
the strength of coupling between the internal motion and 
channels. 

We calculate the average 5'-matrix and distribution 
of complex energy levels in the limit M, iV — > 00, with 
m = M/N < 1 fixed. It turns out that the band struc- 
ture of the hermitean part H does not lead to the essen- 
tial complication when evaluating the one-point average 
characteristics and all calculations can be done in close 
analogy with those of the GOE case Below we 

present rather short description, referring for details to 



3. The calculation of the average iS-matrix is related 
with that of the average Green function 
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which governs the evolution of the intermediate unsta- 
ble chaotic system. It is convenient to use the following 
representation for {Gnm{z)) fH] 
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To carry out averaging in we use the replica method 
[pl|J^]. The generating function Z{z,I) can be repre- 
sented as a multivariable Gaussian integral that makes 
averaging over the Gaussian RBM ensemble trivial. 
The further integration can be performed by means of 
the saddle-point approximation. The saddle-point solu- 
tion turns out to be stable in the replica space and pro- 
portional to Snn'- Therefore, the replicas decouple and, 
as was pointed in H, it is enough to calculate ln(Z(z, J)). 



One has 
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+ ^-(j)n{z + iVV^ - /)„,„(/),„] I , (7) 

where d[(j)] means the product of differentials and above 
equality is valid to irrelevant constant. To make the inte- 
gration over doable, we introduce, following |^,^, new 
variables cr„ = A0^j with the help of ^-functions defined 
by their Fourier representation 



(5(cr„ - A(/)„) = - / d&n exp{--(CT„(T„ - A(/)„d-„ (/)„)} 
TT J 2 

(8) 

After the Gaussian integration over cj) being done, the 
subsequent integration can be performed in the saddle- 
point approximation justified by two large parameters 
N,b ^ 1 . We find after some algebra that the average 
Green function is determined in the following way 



(gnrniz)) = [{z + Xa^.p. + tVV^) 
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by the solution of the saddle-point equations (with re- 
spect to a and a) for the "Lagrangian" C 



^(^o-„(7„ + ln(z + A(T„)) 



--trcln[f + iV"^(z + A6-)-V] , 
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where the diagonal matrix a = diag((Ti, . . . ,(7 at). Note 
that the trace in the last term in (|lO|) runs over the chan- 
nel (M-dimensional) space. 

Similar to the pure RBM case [0, the saddle point 
equations possess the translation invariant (independent 
of n) solution provided that the relation (0) is fulfilled. 
Going along the same line as in |q,P] , we arrive at 



{Qnm{z)) = 



(z - JoX^g{z)){z - JoX^giz) + ijX) 
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Here the notation Jq = J2r ■^(kl) introduced. Due to 
the translation invariance mentioned above the average 
Green function depends on the only combination of tran- 
sition amplitudes. This fact makes the average S'-matrix 
be diagonal and equal to 



.^^ z- JoA^g(z) - ijX 
z — JoX^g{z) + ijX 
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The function g{z) denotes the trace of the average Green 
function which is found to satisfy the cubic equation 



giz){JoX^giz)-z) + l + 
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where the (unique) solution with a negative imaginary 
part has to be chosen 

For a closed system (7 = 0), this cubic equation is re- 
duced to the quadratic one which determines the density 
of the RBM eigenvalues to be given by Wigner's semicir- 
cle law [E2i with the half-radius 



A — \^J Jq , 
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renormalized by the factor a/Jq (this factor reduces to 
unity in the case of the GOE). For an open system, ad- 
ditional rescaling of the coupling constant 
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reduces eqs.(pl|)-(p^) to corresponding ones |^,^ for the 
full matrix case. 

4. The calculation carried out above is valid only for 
the upper half of the complex plane where the Green 
function is analytical, all 5-matrix singularities being lo- 
cated into the lower part (see eq.(||)). The analyticity 
is broken where the density of complex energies differs 
from zero. A special regularization procedure of the pole 
singularities has been proposed in Q to calculate the 
average ^-matrix's pole distribution in the lower half- 
plane of the complex variable z — x + iy. Basing on 
a convenient electrostatic analogy [0,1), the distribution 
of complex levels is considered to be the source of the 
two-dimensional electrostatic field with the potential 

^x,y)^^{\ndet{{z*-n^){z-H) + Sy') . (16) 

The limit <5 — > O"*" should be taken at the very end of cal- 
culations. Applying the two-dimensional Laplacian to $, 
one gets the average density of complex levels (" charges" ) 



AT:p{x,y) = ~A^{x,y) 
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The replica trick can be used again for performing the 
ensemble averaging. The infinitesimal positive S makes 
the matrix in the rhs of eq.([l^) positive definite for any 
z. Therefore, the determinant may be represented as 
the Gaussian integral over a complex TV- vector ■0(1)- 
To make the ensemble averaging possible, the product 
{z* -rO){z -n) is decoupled |J| with the help of 
the complex Habbard-Stratonovich transformation {N- 
vector ■0(2)) ■ One finally arrives at the integral represen- 
tation 



exp{iV$} = / d[0] ( exp{ii/;tA^0'}) 
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where 2N x 2N matrix M. is defined as 
M = 



z* - i5 
i z — H 



introduced 2Af-vectors being ip'^ = ('0(1)-^, ?/'(2)"'") and 
V^'^^ (0(2)^,^(1)^). 



Performing the same steps as in , we find the po- 
tential to be determined by the saddle-point value of 
the "Lagrangian" 



^ ^ Jnni^^ai^n^m^ ^ tr,^ ((7y^(7y^ ) 
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^tra ln(z5 + A(7„) - trQcln[l - iV'^{zs + Xty^Vl] 



Each of N matrices (t„ in (19) has the following structure 



w u 
V w* 



with the real positive u and v whereas (7„ stands for its 
Fourier counterpart. We have also introduced 2N x 2N 
block-diagonal matrices zg obtained from A4 by setting 
there H equal to zero and nm — (<5'n)Q/3<5nm, the 
2M X 2M block-diagonal matrix I is unity in the channel 
subspace and equals to I = diag(l,— 1) in the replica 
subspace, and V is the unit matrix for replica indices 
and corresponds to for others. In the saddle-point we 
have the translation invariant saddle-point equations 



(T = iJaa 

ia .> T ^ 1 • mjXl 

— {zs + iXJqcf) + 1 - I — — — - 

A zs + lAJoa — i^XL 
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which differ from the corresponding equations for the full 
matrix case only in appearing the renormalizing factor 
Jq. Therefore, the explicit solution can be found in our 
case in the same way as it has been done in As a 
result, one concludes that all poles (charges) lie in the 
finite domain of the lower part {y < 0) of the complex 
energy plane defined by the condition x'^ < x^(y) with 
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Inside this region the density of complex energy levels is 
equal to 



47rp(x, J/) 



m 1 — m 



1 
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One can easily see that under rescaling (|14|),(|15D the av- 
erage complex level distribution reduces again to that |9) 
of the full matrix case. 

In conclusion, the band structure of the hcrmitean part 
H results only in the renormalization of the energy scale 
A (|lj) and coupling constant 7 (|l5|) as compared to the 
GOE model In particular, the condition of the "width 
collectivization" [[7|j2^, 7 ~ 1, implies again the natu- 
ral physical condition of the average partial width be- 
ing comparable with the average level spacing |2^]. One 
should expect nontrivial consequences of the band struc- 
ture to appear only while considering the higher correla- 
tion functions. 
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